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ON SOME NON ASYMPTOTIC BOUNDS FOR THE EULER SCHEME 

V. Lemaire 1 and S. Menozzi 2 

Abstract. We obtain non asymptotic bounds for the Monte Carlo algorithm associated to the Euler 
discretization of some diffusion processes. The key tool is the Gaussian concentration satisfied by the 
density of the discretization scheme. This Gaussian concentration is derived from a Gaussian upper 
bound of the density of the scheme and a modification of the so-called "Herbst argument" used to prove 
Logarithmic Sobolev inequalities. We eventually establish a Gaussian lower bound for the density of 
the scheme that emphasizes the concentration is sharp. 
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CO ■ 1.1. Statement of the problem 



1. Introduction 



Let the K d -valued process (X t )t^o satisfy the dynamics 



Xi 



f b(s,X s )ds+ f Ba(s,X s )dW s , (1.1) 
Jo Jo 



. rH . where (Wt)t^o is a d'-dimensional (d! ^ d) standard Brownian motion defined on a filtered probability space 
■ (tt, J-, (J t )oo, P) satisfying the usual assumptions. The matrix B = ( _ ^ d xd ^ is the embedding matrix 

S3 ■ v u (d-d')xd' j 

from R d ' into R d . The coefficients kl+xE^ R d ,a : M+ x R d -> R d ' ® R d ' are assumed to be Lipschitz 
continuous in space, l/2-H61der continuous in time so that there exists a unique strong solution to (1.1). 

Let us fix T > and introduce for (t,x) € [0,T] x R d , Q(t,x) := E[f(T,X% x )], where / is a measurable 
function, bounded in time and with polynomial growth in space. The numerical approximation of Q(t,x) 
appears in many applicative fields. In mathematical finance, Q(t,x) can be related to the price of an option 
when the underlying asset follows the dynamics (1.1). In this framework we consider two important cases: 

(a) If d — d' , Q(t, x) corresponds to the price at time t when X t — x of the vanilla option with maturity T and 
pay-off /. 

(b) If d! = d/2, b(x) = ( J where bi(x) € R d ,62(2;) = ■ • • ,%d')*, Q{t,%) corresponds to the price of 



an Asian option. 



b 2 (x) 
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It is also well known, see e.g. Friedman [Fri75], that Q(t,x) is the Feynman-Kac representation of the solution 
of the parabolic PDE 

d t Q(t, x) + LQ(t, x) = 0, (t, x) e [0, T) x R d , 

Q(T,x) = f(T,x), x e R d , [ ' 

where L stands for the infinitesimal generator of (1.1). Hence, the quantity Q(t,x) can also be related to 
problems of heat diffusion with Cauchy boundary conditions (case (a)) or to kinetic systems (case (b)). 

The natural probabilistic approximation of Q(t,x) consists in considering the Monte Carlo algorithm. This 
approach is particularly relevant compared to deterministic methods if the dimension d is large. To this end 
we introduce some discretization schemes. For case (a) we consider the Euler scheme with time step A := 
T/N, NeW. Set Vi G N, U = iA and for t ^ 0, define <j>(t) = U for U^t< t l+1 . The Euler scheme writes 

X? = x + j\{<t>{s\X^ s) )ds + J*<t{4>(8), Xf (s) )dW s . (1.3) 

For case (b) we define 

X^ = x + J^ bl{ ^^ (1.4) 

where (X^) 1,d := ((A^) 1 , • • • , {X^) d )*. Equation (1.4) defines a completely simulatable scheme with Gauss- 
ian increments. On every time step, the last a" components are the integral of a Gaussian process. 

The weak error for the above problems has been widely investigated in the literature. Under suitable as- 
sumptions on the coefficients b, a and / (namely smoothness) it is shown in Talay and Tubaro [TT90] that 
E D (A) := E x [f(T,X£)] - E x [f(T,X T )] = CA + 0(A 2 ). Bally and Talay [BT96a] then extended this result 
to the case of bounded measurable functions / in a hypoelliptic setting for time homogeneous coefficients b, a. 
Also, still for time homogeneous coefficients, similar expansions have been derived for the difference of the 
densities of the process and the discretization scheme, see Konakov and Mammen [KM02] in case (a), Kon- 
akov et al. [KMM09] in case (b) for a uniformly elliptic diffusion coefficient crcr*, and eventually Bally and 
Talay [BT96b] for a hypoelliptic diffusion and a slight modification of the Euler scheme. The constant C in the 
above development involves the derivatives of Q and therefore depends on /, 6, a, x. 

The expansion of Ed (A) gives a good control on the impact of the discretization procedure of the initial 
diffusion, and also permits to improve the convergence rate using e.g. Richardson-Romberg extrapolation 
(see [TT90]). Anyhow, to have a global sharp control of the numerical procedure it remains to consider the 
quantities 

M 

E MC (M,A) = -]T/(T,(A£r)-E4/(T,A£)] . (1.5) 

i=i 

In the previous quantities M stands for the number of independent samples in the Monte Carlo algorithm 
and ((X t A ) 

t>o)icn mi are independent sample paths. Indeed, the global error associated to the Monte Carlo 
algorithm writes: 

E(M, A) = E D (A) + E M c{M, A), 

where Ed(A) is the discretization error and Emc(M, A) is the pure Monte Carlo error. 

The convergence of £mc(M, A), to when M — > oo is ensured under the above assumptions on / by the 
strong law of large numbers. A speed of convergence can also be derived from the central limit theorem, but 
these results are asymptotic, i.e. they hold for a sufficiently large M. On the other hand, a non asymptotic 
result is provided by the Berry-Esseen Theorem that compares the distribution function of the normalized 
Monte Carlo error to the distribution function of the normal law at order 0(M -1 / 2 ). 

In the current work we arc interested in giving, for Lipschitz continuous in space functions /, non asymp- 
totic error bounds for the quantity Emc(M, A). Similar issues had previously been studied by Malrieu and 
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Talay [MT06]. In that work, the authors investigated the concentration properties of the Euler scheme and 
obtained Logarithmic Sobolev inequalities, that imply Gaussian concentration see e.g. Ledoux [Led99], for 
multi-dimensional Euler schemes with constant diffusion coefficients. Their goal was in some sense different 
than ours since they were mainly interested in ergodic simulations. In that framework we also mention the 
recent work of Joulin and Ollivier for Markov chains [ JO09] . 

Our strategy is here different. We arc interested in the approximation of Q{t, x), t ^ T where T > is fixed. 
It turns out that the log-Sobolev machinery is in some sense too rigid and too ergodic oriented. Also, as far 
as approximation schemes are concerned it seems really difficult to obtain log-Sobolev inequalities in dimension 
greater or equal than two without the constant diffusion assumption, see [MT06]. Anyhow, under suitable 
assumptions on b, a (namely uniform ellipticity of oo* and mild space regularity), the discretization schemes 
(1.3), (1.4) can be shown to have a density admitting a Gaussian upper bound. From this a priori control 
we can modify Hcrbst's argument to obtain an expected Gaussian concentration as well as the tensorization 
property (see [Led99]) that will yield for r > and a Lipschitz continuous in space f,W[\E M c(M,A)\ ^r + S] < 
2 exp(— -^j^r 2 ) for a(T) > independent of M uniformly in A = T/N. Here 8 ^ is a bias term (independent of 
M) depending on the constants appearing in the Gaussian domination (see Theorem 2.1) and on the Wasserstein 
distance between the law of the discretization scheme and the Gaussian upper bound. We also prove that a 
Gaussian lower bound holds true for the density of the scheme. Hence, the Gaussian concentration is sharp, 
i.e. for a function / with suitable non vanishing behavior at infinity, the concentration is at most Gaussian, i.e. 
^[\Emc(M, A)| ^ r — 5\ ^2 exp(— -^yyr 2 ), for r large enough, 5 depending on /, and the Gaussian upper and 
lower bounds, a(T) > independent of M uniformly in A = T/N. 

The paper is organized as follows, we first give our standing assumptions and some notations in Section 1.2. 
We state our main results in Section 2. Section 3 is dedicated to concentration properties and non asymptotic 
Monte Carlo bounds for random variables whose law admits a density dominated by a probability density 
satisfying a log-Sobolev inequality. We prove our main deviations results at the end of that section as well. 
In Section 4 we show how to obtain the previously mentioned Gaussian bounds in the two cases introduced 
above. The main tool for the upper bound is a discrete parametrix representation of Mc Kean-Singer type for 
the density of the scheme, see [MS67] and Konakov and Mammen [KMOO] or [KM02]. The lower bound is then 
derived through suitable chaining arguments adapted to our non Markovian setting. 

1.2. Assumptions and Notations 

We first specify some assumptions on the coefficients. Namely, we assume: 
(UE) The diffusion coefficient is uniformly elliptic. There exists A > s.t. for (t,x,£) G [0,T] x R d x R d ' we 

have Aq 1 |^| 2 < (a(i, x)£, £) < A |£| 2 where a(t,x) := <T<r*(t,x), and |.| stands for the Euclidean norm. 
(SB) The diffusion matrix a is uniformly 77-Holder continuous in space, i] > 0, uniformly in time, and the drift 

b is bounded. That is there exists L > s.t. 



Throughout the paper we assume that (UE), (SB) are in force. 

In the following we will denote by C a generic positive constant that can depend on L$, Ao, 7], d, T. We reserve 
the notation c for constant depending on Lq, \o,rj,d but not on T. In particular the constants c, C are uniform 
w.r.t the discretization parameter A = T/N and eventually the value of both c, C may change from line to line. 

To establish concentration properties, we will work with the class of Lipschitz continuous functions F : R d — > 
R satisfying IVF^ = esssupjVF(x)| < 1 where \VF\ denotes the Euclidean norm of the gradient VF, defined 
almost everywhere, of F. 

Denote now by S^ 1 the unit sphere of M. d . For z £ M d \{0}, ngd-i(z) stands for the uniquely defined 
projection on S 1 . For given po > 0, j3 > 0, we introduce the following growth assumption in space for F in 
the above class of functions: 



sup 

(t,x)£[0,T]xK d 



\b(t,x)\ + 



sup 

te[0,T],(x,y)m 2d , x^y 



a(t,x) - a(t,y)\ 
\x - y\ r ' 
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(G POil g) There exists A C S d 1 such that 

Vy e M d \B( ( o ), 7r sd -i (y) e A, y a := po^-i (l/), ^(l/) - > /% - yo\, 

with A of non empty interior and \A\ ^ e > for d ^ 2 (|.| standing here for the Lebesgue measure of 

and Ac{- 1, 1} for d = 1. In the above equation B(p ) stands for the Euclidean ball of M. d of radius p , and 

(., .) denotes the scalar product in R d . 

Remark 1.1. The above assumption simply means that for \y\ ^ po the graph of F stays above a given 
hyperplane. In particular, for all z e A, F(rz) — > +oo. 

The bounds of the quantities Emc {M, A) will be established for real valued functions / that are uniformly 
Lipschitz continuous in space and measurable bounded in time, such that for a fixed T, F{.) := f(T, .) will be 
Lispchitz continuous satisfying (V-F^ ^ 1. Moreover, for the lower bounds, we will suppose that the above F 
satisfies (G po _p). 

2. Results 

Let us first justify that under the assumptions (UE), (SB), the discretization schemes admit a density. For 
all x G M d , ^ j <j'^N,Ae B(R d ) (where B(R d ) stands for the Borel cr-ficld of R d ) we get 

p [ X v eA \ x u= x ] = / p^{t j ,t j+l ,x,x j+ x)p A {t j+ x,t j+ 2,x j+l ,x j+2 )---p A {tj l -x,t j -,x jl - U x jl ) 

x dxj + idxj + 2 ■ ■ ■ dxji , (2.1) 

where the notation p A (U, Xi, Xi+i), i G [0, N— 1] stands in case (a) for the density at point Xi+i of a Gauss- 
ian random variable with mean Xi + b(ti, Xi)A and non degenerated covariance matrix a(ij, Xi)A, whereas in case 

^ 1 T rl' 1 O I 

x d+1 ' d + x]' d A + hfaxJtfft J 

and non degenerated as well covariance matrix ( a |! i,X ^2 /n ^l^'^us^ )' wnere £ ^ d -> = 

{y 1 ,...,y d ')* and y d '+^ d = (y d ' + \ . . . , y d )* . 

Equation (2.1) therefore guarantees the existence of the density for the discretization schemes. From now 
on, we denote by p A (tj,tj>,x,-) the transition densities between times tj and tj>, ^ j < j' < N, of the 
discretization schemes (1.3), (1.4). Let us denote by P x (resp. Pt. x , ^ j < N) the conditional probability 
given {X^ = x} (resp. {X^ = x}), so that in particular F X [X£ E A] = J A p A (0,T,x,x')dx' . We have the 
following Gaussian estimates for the densities of the schemes. 

Theorem 2.1 ("Aronson" Gaussian estimates for the discrete Euler scheme). Assume (UE), (SB). There 
exist constants c > 0, C ^ 1, s.t. for every < j < j' ^ iV: 

C~ 1 p c -i(tj> —tj,x,x') < p A (tj,tj> ,x,x') ^ Cp c (tj> —tj,x,x'), (2.2) 

/ \ d/2, i / 1 2 

where for all ^ s < t < T , in case (a), p c (t — s,x,x') := ( 2 -K{t-s) ) ex P(~ c 't(t-s) ) ant ^ * n case W 

„ / VSc Y /2 ( \\{x'Y' d ' - x^ d ' \ 2 Ux'f +1 - d - x d ' +1 ' d - * 1 - d '+(*') 1 - d ' (t - s)\ 2 

*<t->*>^ ■= {^-s? ) cx pH \ { t- s) +3 — (i^r^ — 

Note that p c enjoys the semigroup property, i.e. V0 < s < t, J Rd p c (t — s, x, u)p c {s, u, x')du — p c (t, x, x') (see 
Kolmogorov [Kol34] or [KMM09] for case (b)). 
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Remark 2.1. The above upper bound can be found in [KM02] in the case of time homogeneous Lipschitz 
continuous coefficients. Both bounds can be derived (for time dependent coefficients) from the work of Gobet 
and Labart [GL08J under stronger smoothness assumptions. Here, our framework is the one of the "standard" 
PDE assumptions to derive Aronson's estimates for the fundamental solution of non degenerated non- divergence 
form second order operators, see e.g. Sheu [She91] or [DM09]. In particular no regularity in time is needed. 

Our second result is the Gaussian concentration of the Monte Carlo error Emc(M, A) defined in (1.5) for a 
fixed M uniformly in A = T/N, N > 1. 

Theorem 2.2 (Gaussian concentration). Assume (UE), (SB). For the constants c and C of Theorem 2.1, we 
have for every A = T/N, N ^ 1, and every Lipschitz continuous function in space and measurable bounded in 
time f : R d -> K satisfying |V/(T, .J^ ^ 1 in (1.5), 



Vr>0, VM>1, F x [| E MC (M, A)| ^ r + 5 c , a ( T )] < 2e" 



(2.3) 



with 



a(T) 




in case (a), 
^ in case (b), 



(2.4) 



and S c , a (T) = 2 v /a(T) logC. 

Moreover, if F(.) := f(T, .) ^ satisfies for a given p > and (3 > 0, the growth assumption (G Po jj), 



Vr > 0, VM > 1, P x [|£ M c(M,A)| > r - <J c ,c, r ,/] > 2exp 



/ M 


r 


I a(T) 





(2.5) 



where 5 c ,c,t,/ = (1 + V2)y/a(T) logC + 7 c -i iT (F) + p /3 - £, 7 c -i, T (da; / ) = p c -i (T, x, x')dx' , and F := 
inf seS d-i F(spo). The constant a(T) _1 appearing in (2.5) writes in case (a) 



1 



a(T) 



A + x := 



2T 



2T 



c 

K(d,A) 



+ 



log(- 



" y arccos(e- 1 /2)K((i, J 4) / + 

w/iere /or all d e N* , Ac S"^ 1 appearing in (G po ^), 



for d even, 

for d odd, 9 e (1, +oo), 



K(d, A) = 



|A|(d/2-l)! 



7rV2 



d even, 

d odd. 



(2.6) 



In case (b), d is even and 



1 t c- 1 / 3 

= A + y := H 

a{T) X 2T \ T 2 







2^4 




1 + 






M iog ( 



(f) d/2 [^ 2 + 3(l + v/l 



3 



A"(d,A) 



T' 1 



)] d / 2 c N 



From Theorem 2.1 and our current assumptions on /, we can deduce from the central limit theorem that 

M 1 /2 £ ; MC (M,A) AA(0,a 2 (/,A)), a 2 (/,A) := E x [f(Xfi) 2 } - E x [f(Xfi)} 2 . From this asymptotic regime, 

we thus derive that for large M the typical deviation rate r (i.e. the size of the confidence interval) in (2.3) has 
order ca{f, A)M -1 / 2 where for a given threshold a e (0, 1), c := c(a) can be deduced from the inverse of the 
Gaussian distribution function. In other words, r is typically small for large M. On the other hand, we have 
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a systematic bias <5c,a(T)> independently of M. In whole generality, this bias is inherent to the concentration 
arguments used to derive the above bounds, see Section 3, and cannot be avoided. Hence, those bounds turn 
out to be particularly relevant to derive non asymptotic confidence intervals when r and &c.a{T) have the same 
order. In particular, the parameter M is not meant to go to infinity. This kind of result can be useful if for 
instance it is particularly heavy to simulate the underlying Euler scheme and that only a relatively small number 
M of samples is reasonably allowed. On the other hand, the smaller T is the bigger M can be. Precisely, one 
can prove that the constant C of Theorem 2.1 is bounded by cexp(L ^) (see Section 4). Hence from (2.4), we 
have 5 c>a (T) = 0(T) for T small. 

Remark 2.2. For the lower bound, the "expected" value for <5(T) -1 would be A corresponding to the largest 
eigenvalue of one half the inverse of the covariance matrix of the random variable with density p c -i(T, x, .) 
appearing in the lower bound of Theorem 2.1. There are two corrections with respect to this intuitive approach. 
First, there is in case (a) an additional multiplicative term, 9 > 1 (that can be optimized) when d is odd. This 
correction is unavoidable for d = 1, anyhow for odd d > 1, it can be avoided up to an additional additive factor 
like the above \ (see the proof of Proposition 3.3 for details). We kept this presentation to be homogeneous for 
all odd dimensions. 

Also, an additive correction (or penalty) factor \ appears. It is mainly due to our growth assumption 
(G POt p). Observe anyhow that, for given T > 0,C ^ l,e > s.t. \A\ > e, if the dimension d is large enough, by 
definition of K(d,A), we have \ = 0. Still, for d = 1 (which can only occur in case (a)) we cannot avoid the 
correction factor \- 

Remark 2.3. Let us also specify that in the above definition of x, Pa is not meant to go to zero, even though 
some useful functions like |.| satisfy (G Po ,i) with any po > 0. Actually, the bound is particularly relevant in 
'large regimes", that is when r/(3 is not assumed to be small. Also, we could replace in the above definition of x, 
Po by R > as soon as r/f) ^ R. In particular, if F satisfies (G po / g) ; for i? ^ po it also satisfies (Gr^). We 
gave the statement with po in order to be uniform w.r.t. the threshold po appearing in the growth assumption of 
F but the correction term can be improved in function of the deviation factor r/(3. 

Remark 2.4. Note that under (UE), (SB) 7 in case (a), the martingale problem in the sense of Stroock 
and Varadhan is well posed for equation (1.1), see Theorem 7.2.1 in [SV79]. Also, from Theorem 2.1 and 
the estimates of Section 4, one can deduce that the unique weak solution of the martingale problem has a 
smooth density that satisfies Aronson like bounds. Furthermore, a careful reading of [KM02] emphasizes that 
the discretization error analysis carried therein can be extended to our current framework, i.e. we only need 
boundedness of the drift and uniform spatial Holder continuity of the (non-degenerated) diffusion coefficient to 
control Ed(A). Hence, the above concentration result gives that in case (a), one can control the global error 
E(M,A) := Ed (A) + E M c(M, A). The well-posedness of the martingale problem in case (b) remains to our 
best knowledge an open question and will concern further research. 

Remark 2.5. In case (b), the concentration regime in the above bounds highly depends on T. Since the two 
components do not have the same scale we have that, in short time, the concentration regime is the one of the 
non degenerated component in the upper bound (resp. of the degenerated component in the lower bound). For 
large T , it is the contrary. 

We now consider an important case for applications in case (b). Namely, in kinetic models (resp. in financial 
mathematics) it is often useful to evaluate the expectation of functions that involve the difference of the first 
component and its normalized average (which corresponds to a time normalization of the second component). 
This allows to compare the velocity (resp. the price) at a given time T and the averaged velocity (resp. averaged 
price) on the associated time interval. Obviously, the normalization is made so that the two components have 
time-homogeneous scales. We have the following result. 



Id'xd' Od'xd' 
Od'xd' r _1 Id'xd 

and g is a Lipschitz continuous function in space and measurable bounded in time satisfying | Vt/(T, .) loo ^ 1 



Corollary 2.1. In case {b), if f in (1.5) writes f(T,x) = g{T,T^}x) where T T X 
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then we have for every A — T/N, N > 1, 

VM >1, P x [|£7mc(M,A)| >r + 5 CMT) ] ^2e 



with a{T) = (4 - >/l3)# £c,a(T) = V"^) logC. 



A lower bound could be derived similarly to Theorem 2.2. 

The proof of Theorems 2.1 and 2.2 (as well as Corollary 2.1) are respectively postponed to Sections 4.2 and 



3.3. 



3. Gaussian concentration and non asymptotic Monte Carlo bounds 



3.1. Gaussian concentration - Upper bound 

We recall that a probability measure 7 on R d satisfies a logarithmic Sobolev inequality with constant a > 
if for all / e H 1 ^) := {g G L 2 (d 7 ) : / |V. 9 | 2 d7 < +00} such that / > 0, one has 



Ent 7 (/ 2 ) < a J |V/| 2 d 7 , 



(CS2 a ) 



where Ent 7 (</>) = / 01og(0)d7 — ( f </>d 7 ) log ( f 0d 7 ) denotes the entropy of the measure 7. In particular, we 
have the following result (sec [Lcd99] Section 2.2 eq. (2.17)). 

Proposition 3.1. Let V be a C 2 convex function on R d with HessV^ ^ Aldxd; A > and such that e _1 



is integrable with respect to the Lebesgue measure. Let 7 (dx) 



vOOda; be a probability measure (Gibbs 



measure). Then 7 satisfies a logarithmic Sobolev inequality with constant a = j. 

Throughout this section we consider a probability measure fi with density m with respect to the Lebesgue 
measure Xk on (here we have in mind K = d or K = Md, M being the number of Monte Carlo paths). 
We assume that /i is dominated by a probability measure 7 in the following sense 



7 (dx) = q(x)dx satisfies (CSl a ) and 3k ^ 1, Vi e R K , m(x) ^ nq(x). 

Proposition 3.2. Assume that fj, and 7 satisfy (H K a ). Then for all Lipschitz continuous function F : 
s.t. IVFU < I, 



Vr > 0, 



F(Y)-n(F)^r + W 1 (ji,'r) 



C ne 



where Wi(fi, 7) = sup \fi(F) — j(F)\ (Wasserstein distance W\ between \i and 7 ) . 

IVFI^l 

Proof. By the Markov inequality, one has for every A > 0, 



F(Y)-^(F)^r + W 1 (p, 1 ) 



< e -A(p(F)+r+W 1 ( M , 7 )) E 



,\F(Y) 



(3.1) 



and by (H K , a ), E M [e AF ( y )] ^ KE 7 [e AF ^ y - 1 ] . Since 7 satisfies a logarithmic Sobolev inequality with constant 
a > 0, the Herbst argument (see e.g. Ledoux [Lcd99] section 2.3) gives 



E- 



,XF(Y) 



< e A7(F) + f A 2 



so that E (U [e AF ( r )] < Ke AMF)+f a^+a^f)-^)) and 



E, 



s AF(y) 



< K e v(F)+ ^ A2+Am( ' i ' 7) 



(3.2) 
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since owing to the definition of W\ one has Wi(/j,, 7) > 7(F) — /i(F). Plugging the above control (3.2) into (3.1) 
yields 

P^F(F) - (1(F) ^ r + Wi(m, 7)] ^ Ke - Ar +? A2 . 
An optimization on A gives the result. □ 
Lemma 3.1. Assume that \i with density m and 7 with density q satisfy the domination condition 

3k > 1, Va; € M d , m(x) ^ nq(x) 

and that there exist (a, ft, ft) € (R+) 3 such that for all Lipschitz continuous function F satisfying (VF^ ^ 1 
and for all A > ; E 



,\F(Y) 



<J e *l(F) + f A 2 +/3iA+/3 2 _ Tftgrl wg ^(^^ <; ft + (ft + log( K )). 



Proof. Recall first that for a non- negative function /, we have the following variational formulation of the 
entropy: 

Ent 7 (/) = sup {E 7 [//i] ; E 7 [e /l ] 1} . (3.3) 
W.l.o.g. we consider F such that fj,(F) ^ 7(F). Let A > and h := AF - A 7 (F) - f A 2 - ft A - ft so that 
E 7 [e' 1 ] 1 and 



E„[A] = E n 



m 



A( M (F)- 7 (F))--A 2 -ftA- 



We then have 



A1 (F)-7(F) = |A + ft + i(^ft+E 7 

(3.3) Q, 1 / 

< T A + ft + T ft 



— h 



Ent 7 I — 



An optimization in A yields 

MF) - 7(^) < ft 
Now using the domination condition, one has Ent 7 ( ™ 




a ft + Ent 7 — 



(3.4) 



/ y log ( y J d7 ^ log(/c) and the results follows. □ 



Remark 3.1. 7Vo£e that if 7 satisfies an (LSI a ) we have ft = ft = and i/ie result (3.4) of Lemma 3.1 

reads Wi(fi, 7) ^ y^aEnt-y ^) ^ \/alog(K). For similar controls concerning the W2 Wasserstein distance see 

Theorem 1 of Otto and Villani [OV00] or Bobkov et al. [BGL01]. 

Using the tensorization property of the logarithmic Sobolev inequality we derive the following Corollary. 
Note that the term S K ^ a can be seen as a penalty term due on the one hand to the transport between \i and 
7, and on the other hand to the explosion of the domination constant n M between /z® M and 7® M when M 
tends to infinity. We emphasize that the bias 5 K , a is independent of M. Hence, the result below is especially 
relevant when r and 5 K , a have the same order. In particular, the non-asymptotic confidence interval given by 
(3.5) cannot be compared to the asymptotic confidence interval deriving from the central limit theorem whose 
size has order 0(M -1 / 2 ). 

Corollary 3.1. Let Y 1 , . . . ,Y M be i.i.d. M. d -valued random variables with law \x. Assume there exist a > 
0, k > 1 and 7 such that {H Kt<x ) holds on R d . Then, for all Lipschitz continuous function f : M. d — )■ M satisfying 

l V /loo ^ 1 > We haVe 



Vr > 0, M > 1, 



1 M 



k=i 



< 2e 



-M — 



(3.5) 



NON ASYMPTOTIC BOUNDS FOR THE EULER SCHEME 



9 



with 5 K . a — 2\fa log(/c) ^ 0. 

Proof. Let r > and M > 1. Clearly, changing / into — /, it suffices to prove that 



1 M 



fe=i 



^ e 



-M — 



By tensorization, the measure ~f® M satisfies an (£SI a ) with the same constant a as 7, and then the probabilities 
/z® M and 7® M satisfy (K« a) on M. K ,K = Md. In this case, Lemma 3.1 gives y/MS K , a ^ Wi(/x® M ,7® M ) + 
y/ Ma log(«;) and then 



fe=i 



, M _ 



Applying Proposition 3.2 with the measures [i® M and 7 (glM , the function F(x±, . . . , xm) = "^7 2~^fc=i 
(which satisfies (VF^ < I) and f = \/~M(r + ^/alog(K)) we obtain 



M 



M 



k=l 



< K M e~ M 



(r+Volog(«)) ! 



and we easily conclude. 



□ 



Remark 3.2. Note that to obtain the non- asymptotic bounds of the Monte Carlo procedure (3.5), we successively 
used the concentration properties of the reference measure 7, the control of the distance W\ (/1, 7) given by 
the variational formulation of the entropy (see Lemma 3.1) and the tensorization property of the functional 
inequality satisfied by 7. The same arguments can therefore be applied to a reference measure 7 satisfying a 
Poincare inequality. 

3.2. Gaussian concentration - Lower bound 

Concerning the previous deviation rate of Proposition 3.2, a natural question consists in understanding 
whether it is sharp or not. Namely, for a given function / satisfying suitable growth conditions at infinity, 
otherwise we cannot see the asymptotic growth, do we have a lower bound of the same order, i.e. with Gaussian 
decay at infinity? The next proposition gives a positive answer to that question. 

Proposition 3.3. Let f : R d — > M + be a Lipschitz continuous function satisfying IV/I^ ^ 1 and assumption 

{Gpo,p) f or 9 iven Po,P > 0. 

For a C 2 function V on M. d such that e~ v is integrable with respect to Ad and s.t. 3X ^ 1, Aldxd ^ Hess(y) 
0, ie^7(da;) = e~ v ( x > Z _1 dx be the associated Gibbs probability measure. We assume that 3k ^ 1 s.t. for \x\ ^ p 
the measures n(dx) = m(x)dx and j(dx) satisfy 

m{x) > K- x e~ v ^Z- x . 

Let A •= - + sup °es d - 1 I v ( s p°)I _i_ su P a£S ''-i \w(spoV 



+ 



10 



V. LEMAIRE AND S. MENOZZI 



We have 



Vr>0, P M [/(F) -/*(/) >r-(W 1 ( f i, 1 )+S(f, 1 ))} > { 



K(d,A) / T 

arccos(6>~ 1/2 )_ff(tlA) 
ZA d / 2 K 



cxp -OA 



d even, 



Vp 



, V0 > 1, d odd, 



with S(f, 7) = -f(f) + (3po — /, f := ini seS d-i f(sp ), and K(d,A) defined in (2.6) where A C S" i_1 appears in 

( G Po,/3)- 

Proof. Set £ := {y € 4 x [p , 00)}. Here we use the convention that for d = 1, A x [p , +00) C (—00, — po] U 
[poi +00). Write now 



P^[/(F) - ^ r - (Wi(m,7) + «(/,7))] > WOO - M/) > r - (Wi(/x, 7 ) + <K/, 7 )),£] 

^ «- x P 7 [/(F) >r-pp +l, E] := k" 1 ?. 

Denoting Y = p TT S d-i(Y), we have 



(3.6) 



V > P 7 [/(F ) + (/(F) - /(F )) ^r-Ppo+l, E] 



( G PO.f3) 



> P 7 [/3|F-F | >r-0po + f-f{Y o ),E] > 



\Y\> r -^ + \Y \,E 



|F| > - Vpo,7r s «i-i(y) G A 



(3.7) 



Write 



r r+cc 

V > <r(ds) dp/- 1 exp(-F(sp))Z- 1 , 

JA JpnV^ 



'Po 

where cr(ds) stands for the Lebesgue measure of S . Now, Hess(F) ^ Aldxd yields Vp ^ po V ^, \V(sp)\/p 2 < 

and therefore 



^ A -|- SU P 3 gS d -l \ V ( S P0)\ _|_ SUP a g S d-l |W(sp )| 



V > \A 



Po 



7 PoV ? 
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JpnV^ Z(2A) / J( po v^)(2A)i/ 

^(p V^)(2A) 1 ^. 



Z(2A)<V 2 

We now have the following explicit expression: 



lAl 1 ~ dpp^cM-^) 



(3.8) 



Vx > 0, Q d (x) 



M(d,x) 



exp(-y)M(d,x), 



y^.r' TT 2j, d even, 



i=0 j=i+l 



2 x 2l+1 [] (2j + 1) + Yl (2j + 1) cxp(-) / cxp(-^)dp, d odd, 

i=0 J=t 3=0 1 Jx 1 



with the convention that J2i2 = 0, Vfc € N, IlLfe J = L 
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Observe now that /~ exp(-p 2 /2)dp = (2tt) 1 / 2 P[JV(0, 1) > x] > (2 7 r) 1 / 2 P[r e K,\Y\ > a;/cos(<9)] := 
(2 7 r) 1 / 2 Q(a;), where F ~ A/"(0 2 xi, I2X2) is a standard bidimensional Gaussian vector and JC := {z € R 2 , (z, ei) > 



cos(0)|z|}, 9e (0,f), ei = (1,0). Since = f exp(- 



2cos 2 (#) 



), we derive that 



i{x) > 



2<*/2-i(d/2- l)!exp(-^), dcven, 



S^n^(i-|)exp(-^M dd, 



2c OS 2(e) • 



which plugged into (3.8) yields: 



y ( d ^) exD [ -A 

ZA d / 2 ^ 1 



0if(d.A) 

^T^exp 



A 



COS 2 (9) 



d even, 
2^ 



d odd. 



□ 



Corollary 3.2. Under the assumptions of Proposition 3.3, let Y 1 , • • • , Y M be i.i.d. M. d -valued random variables 
with law jj,. We have Vr > 0, VM ^ 1, 



2 x < 



exp f -M(A + x) 
exp ( -M(9A + x) 



^E/C^-M/O' 1 )] 
fc=i 



>r-(Wi(/i, 7 )+<f(/,7)) 



,A-(d,A)arcco^T72); + ' ^ ^ (1 , +Oo) , d orfrf, 

wzi/i if(d, A) defined in (2.6). 

Proof. We only consider <i even. By independence of the {{Y) k )keli.M\ > exploiting rifcli{/(^ fe ) — [/C^ 1 )] ^ 
r-(l^ 1 ( M ,7)+<5(/,7))}c{^Ef = i/(^ fe )-E4/(y 1 )] ^-(^ )7 )+<S(/,7))},wehave 



^ E /( yfe ) - ^ > r - 7) + *(/, 7)) 

fe=i 

For x = ^ log ( f^) ) ' we thus obtain 



( K(d, A) 
\ ZA d / 2 n 



M 



exp -MA 



Vp 



1 M 

- J2 f(Y k ) E M [/(F 1 )] > r - (Wi (/x, 7 ) + <*(/, 7)) 



fe=i 



^ exp — M(A + x) 



which completes the proof. 



3.3. Proofs of Theorem 2.2 and Corollary 2.1 



□ . 



- Theorem 2.2 - Upper bound (2.3). 

In case (a), the Gaussian probability 7 Ci t with density p c {T, x, .) dchned in Theorem 2.1 satisfies a logarithmic 
Sobolev inequality with constant a(T) — ^L. The result then follows from Theorem 2.1 and Corollary 3.1. 
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In case (6), 7 c ,T(dx') = p c (T,x,x')dx' = Z- x e~ VT ^dx' where 



(~r\l,d' „l.d'\2 \( T >\d' + l,d „d'+l,d x +(x ) T |2 

y T ,,(^) = c( l( } _ 1 + 3 K j 2 11 



AT 



T 3 



(3.9) 



The Hessian matrix of Vts satisfies 



Vx' G M d , Hess^x') = f &' xd ' i Id ' xd ' ) > AI dxd , 

\ -J^--ld'xd' T^T-ld'xd' / 



with A = f + |§ 



1 — -v/l + ^ + ^l > 0. By Proposition 3.1, the probability 7 Ci t satisfies a logarithmic 



2T 



-. We still conclude by Theorem 2.1 and 



Sobolcv inequality with constant a(T) = 7 ? — — 

Corollary 3.1. 

- Theorem 2.2 - Lower bound (2.5). 

With the notation p c -i(t — s,x,x') = Z^ 1 e~ Vt - £ - x ^ x \ the Hessian of the potential Vt, x satisfies Vx' € M d , 

HessVr,:i;(a;') ^ AIdxd where A = in case (a) and A = + ^A ^1 + yA + ^ + in case (6). 

Set 7 c -i j r(da; / ) = p c -i(T, x, x')dx' and fiT(dx') = p A (0, T, x, x')dx' . Since /ir and 7 c ,t satisfy (H Kt a) with 
k = C and a = a(T) defined in (2.4), the probability [it satisfies (3.2), and Lemma 3.1 yields W\{p,T, 7c,t) ^ 
^/ a(T) log(C). Now, 7 c -i t an d 7c, t satisfy {T-L K ,a) with k = C 2 and a = a(T). We therefore get from Lemma 
3.1, V Ki(7c-i,t,7c ,t) < V / 2«( T ) 1 °g( C ')- _ Hence, Wi(7 c -i, r , A*r) < Wi(/ir, 7c,t) + Wi(7c-i,r, 7c,r) < (1 + 
\/2) ^/a^T) l°g(C)- Now , by definition of 6 c ,c,T,f we have S c> c,T,f > Wi(~f c -i, T , fj, T ) + £>(/> 7c-\t), {S(f, 7c- 1 ,t) 
introduced in Proposition 3.3) and Corollary 3.2 yields 



M 



M ^ 

fe=i 



/(T, (4)*) - [/(T, X A )] ^ r - -V, •. / . / 



^ exp 



M 
~W) 



(3 



where observing that for our Gaussian bounds A = -|, and 



X, X = jz log / 2 

Pa 



- d / 2 1 2 

AT(<Z,-A) 



for d even, 



^ H+x,x = A 1 °g( ^=^»y ) + *»dodd,*>i, 



and defined in (2.6). 

Observe now that in case (a), the normalization factor Z = Z(T,d) associated to p c -i(T, x, .) writes Z 
(2ircT) d / 2 . Hence, recalling that A = (cT) _1 , we obtain in this case 



^ lo g(flS)) + ford even, 
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In case (b), wc have Z = (2^c) d / 2 T d , A = ^ (\ + £ 1 + ^1 + ^ + ^ ^ so that 2- d / 2 ZA d / 2 = (f ) d/2 [T 2 + 
3(1 + v /l + ^ + ^)] d / 2 . Eventually, since in case (b) we always have d even, the correction writes 

1 / (f) d/2 [T 2 + 3(l + v /l + f + ^ C S 
X Pi g \ K{d,A) 

This completes the proof. □ 
- Proof of Corollary 2.1. 

Note that the random variable Yfi = Y T X X£ admits the density p$(T, y, y') = T d 'p A (0, T, T T y, T T y') with 
respect to Xd(dy'). By Theorem 2.1 this density is dominated by (ZT d )- 1 e ~ VT - T Tv( T Ty ) w here Vt, x is defined 
in (3.9). The Hessian of y' i-> Vtj tV (TtZ/') satisfies 

\/y' G M d , HessV+, TTy (T T y') = f J 3 j d ' xd ' 1t* d '* d ' ) > AI dxd , 



□ 



-Id'xd' ^Id'xd 

with A = y(4 — \/l3). We still conclude by Proposition 3.1, and Corollary 3.1. 



4. Derivation of the Gaussian bounds for the discretization schemes 

4.1. Parametrix representation of the densities 

We first derive a parametrix representation of the densities of the schemes. The key idea is to express this 
density in terms of iterated convolutions of the density of a scheme with frozen coefficients, that therefore 
admits a Gaussian density, and a suitable kernel, that has an integrable singularity. These representations have 
previously been obtained in Konakov and Mammen [KMOO] and Konakov et al. [KMM09] . 

We first need to introduce some objects and notations. Let us begin with the "frozen" inhomogeneous scheme. 
For fixed x,x' e R d , < j < f < N, we define (X*) i£[ jl} (= {x£ X ') ieVJ , } ) by 

X* =x, Vie [j,f), X£ +1 = X* + b(U,x')A + a(ti,a/)(W ti+1 - W u ) (4.1) 

for case (a). Note that in the above definition the coefficients of the process are frozen at x', but we omit this 
dependence for notational convenience. In case (b) we define (AQ A ) jg (= (X^' x by X A = x, and 

v/elj./). 

t.+ l I , n /+ ' / Vd'xl 




: = ^ 1- I / ! * +1 ( ^A )M ' ds J +Sa ^ ( ^> ) (tf-U)) (W ti+1 -W ti ). (4.2) 

That is, in case (b) the frozen process also depends on j' through an additional term in the diffusion coefficient. 
This correction term is needed, in order to have good continuity properties w.r.t. the underlying metric associ- 
ated to p c when performing differences of the form a(tj,x) — a(tj,x' — ^ ^ (*?' — *»))) sec * ne definition 

(4.7) and Sections 4.2 and 4.3 for details. 

From now on, p A (tj , ij' , x, •) and p A,< j' ,a: (ij , ty , x, •) denote the transition densities between times tj and ty 
of the discretization schemes (1.3), (1.4) and the "frozen" schemes (4.1), (4.2) respectively. 
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Let us introduce a discrete "analogue" to the inhomogeneous infinitesimal generators of the continuous objects 
from which we derive the kernel of the discrete parametrix representation. For a sufficiently smooth function 
ip : R d -» R and fixed x' G R d , f G (0,iV], define the family of operators (£ A )je[o,.j') and {L A )jelo,j') (= 

rA , n E[^(X A +A )|X A =x]-#r) E[^X A +A )|X A = *]-</>(*) 
= " ^ , and L t .V(a:) = 1 ^ • 

Using the notation p A (tj,tj>,x,x') = p^V' 1 (tj,tji ,x,x'), we now define the discrete kernel H A by 

H*(t j ,t j ,,x,x?) = (Lf;-Lf^^(tj+A,t jt ,x,xf), 0^j<j'^N. (4.3) 

Note carefully that the fixed variable x' appears here twice: as the final point where we consider the density 
and as freezing point in the previous schemes (4.1), (4.2). Note also that if / = j + 1 i.e. ty = tj + A, the 
transition probability p^'*j'' x (tj+i,tj+i, ., x 1 ) is the Dirac measure 5 X > so that 

H A ( tj ,t j+1 ,x,x') = A- 1 (e[5 x ,(X a +i )\X a = x] -E[6 X ,(X A +1 )\X A = x}) , 
= A" 1 (p A {tj,t j+1 ,x,x') -p A ^" x \t j ,t j+1 ,x,x') S j . 

From the previous definition (4.3), for all < j < j' ^ N, 

H A (tj,t r ,x,x') = A -1 / [p A -p^'V'*'] (tj-.^+i.^.u^'V^'^+^tj-^ti^Odu. 

Analogously to Lemma 3.6 in [KMOO] we obtain the following result. 

Proposition 4.1 (Parametrix for the density of the Euler scheme). 
Assume (UE), (SB) are in force. Then, for < t j < tj' ^ T, 

j'-j 

p A (tj,tj,,x,x') = J2 (p A ®a # A ' W ) (tj,tj,,x,x'), (4.4) 

where the discrete time convolution type operator (g) A is defined by 

j'-j-i 

(0Oa f){tj,tj,,x,x') = V A / g(tj,tj +k ,x,u)f(tj + k,tj>,u,x')du, 

k=o jRd 

where g ®a H A '^ = g and for all r ^ 1, H A '^ = H A ®a H A '( r ~ 1 ^ denotes the r-fold discrete convolution of 
the kernel H A . W.r.t. the above definition, we use the convention that p A ®a H A '( r \tj,tj,x,x') — 0,r ^ 1. 

4.2. Proof of the Gaussian estimates of Theorem 2.1 

The key argument for the proof is given in the following lemma whose proof is postponed to Section 4.3. 
Lemma 4.1. There exists c> 0, C ^ 1 s.t. for all < j < f TV, /or aZZ r G [0,./ - j], \/(x, x') G R d , 

\p A ® A ff^fe, t f ,x,x')\ < C r+1 (V - t/"/ 2 'f[B (l + ^=2^, f) Pc(V - tj,x,x'). (4.5) 

Z— 1 ^ ~ ' 
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In the above equation B{m 1 n) := L s m_1 (l — s)™ _1 ds stands for the (3 function. 

The upper bound in (2.2) then follows from Proposition 4.1 and the asymptotics of the /? function. It is also 
useful to achieve the first step of the lower bound. 

Proof of the lower bound. We provide in this section the global lower bound in short time. W.l.o.g. we 
assume that T ^ 1. This allows to substitute the constant C appearing in (4.5) by a constant c < cexpd&loo) 
uniformly for ty — tj ^ T. From the upper bound, we derive the lower bound in short time, on the compact 
sets of the underlying metric, see (4.7) below. This gives the diagonal decay. To get the whole bound in short 
time it remains to obtain the "off-diagonal" bound. To this end a chaining argument is needed. In case (a) it is 
quite standard in the Markovian framework, see Chapter VII of Bass [Bas97] or Kusuoka and Stroock [KS87] . 
In case (b), the chaining in the appendix of [DM09] can be adapted to our discrete framework. We adapt below 
these arguments to our non Markovian setting for the sake of completeness. 

Eventually, to derive the lower bound for an arbitrary fixed T > it suffices to use the bound in short time 
and the semigroup property of p c -i. Naturally, the biggest is T, the worse is the constant in the global lower 
bound. 

From Proposition 4.1 we have 

j'-j 

p A (tj,t J ,,x,x / ) ^ p A (tj,tf,x,x') - \p A <E>a H A ' {r \tj,t r ,x,x')\ 

^ CQ 1 p c -i(t r - tj,x,x') - c (ij' - tj) v/2 p c (tj> -tj,x,x'), (4.6) 

exploiting p A {tj, tj>, x, x 1 ) ^ c " 1 p c -i (£,■/ —tj,x,x') (cf. Lemma 3.1 of [KMM09] in case (b)) and (4.5) (replacing 
C by Co) for the last inequality. Equation (4.6) provides a lower bound on compact sets provided that T is small 
enough. Precisely, denoting 



l -: m case (a), 

(4.7) 



CO 



l,d' „l,d' 



2(V-*j) 



+ 6 



(x') d ' +1 - d -x d ' +1 - d - xl ' d ' + l x,)1 ' d ' (ht-tj) 



(tj'-tj) 



in case (b), 



2/r, 



we have that, for a given R ^ 1/2, if d 2 t _ t . (x, x') < 2R and (tj> — tj) < T < (c 2 exp(— c 1 R )/2) 

1 c" 1 
p A {t j ,t j ,,x,x') > 7 -__( c - 1 exp(- C - 1 i?o)-cor"/ 2 ) ^ °_ cxp(-c- 1 R ) 

(tj> tj) Z(tj> tj) 

where the parameter S is the intrinsic scale of the scheme. In case (a) S = d/2, in case (b) S = d. Hence, up 
to a modification of Cg 1 we have that 

-l 

3c > 1, VO < j < f ^ N, V(x,x') e (M. d ) 2 , d 2 t .(x,x') < 2R , P A (t J ,t r ,x,x r ) > % (4.8) 

v-j' t j! 

In particular 3c > 0, c > 1,V0 < j < j' < N, M{x,x') <= (M d ) 2 , <% t} {x,x') 2i? , p A (^-, , x, a;') ^ 

Cq 1 p c -i(t jl - tj,X,x'). 

Chaining in case (a). Let us introduce: VO s < t sC T, {x,x' ,y) G (R d ) 3 , 2J A ' v (s, i, x, x')Ax' := ¥[Xf e 
dx'|X A = i,l^ s j = y]. Equation (4.8) provides a lower bound for the density of the scheme when s,t 
correspond to discretization times. For the chaining the first step consists in extending this result to arbitrary 
times < s < t < T. Precisely, if d 2 _ s (x, x') < Rq/12 we prove that 

3c > 1, VO < s < t ^ T, Vy, p A ' v (s, t, x, x') ^ c^ \t - s)- d/2 . (4.9) 
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If (p(t) — 4>(s), the above density is Gaussian and (4.9) holds. If 4>{t) = ((f>(s) + A), equation (4.9) directly 
follows from a convolution argument between two Gaussian random variables. Note anyhow carefully that the 
"crude" convolution argument cannot be iterated L times for an arbitrary large L. Indeed, in that case the 
constants would have a geometric decay. Thus, for <j>{t) — (</>(s) + A) > A we write 



p A ' y (s,t,x,x') = / p A ' v (s,(t)(s) + A,x,x 1 )p A ((f>(s) + A, (j)(t),xi,x 2 )p A (4>(t),t,X2,x')dx 1 dx 2 

> I p A ' v {s,<j){s) + A, x, x l )p A {(j)(s) + A,</>(i),xi,x 2 )p A (</>(i),t,x 2 ,x')dxidx 2 

J B R (s,t,x,x') 



Bfi(s,t,x,x') 

(4.10) 

where B R (s,t,x,x') := {x\ G R d : d^,, +A _ s (x, xi) s? R} x {x 2 G M. d : d 2 _ 4> , t Jx 2 ,x') < R} for R > to be 
specified later on. Now, for (x\,X2) G Bn(s,t,x,x'), 

2 \xi -x 2 \ 2 2jgi - x| 2 + 4|x - xf + 4|x 2 - x'| 2 



0(t) -(>(*)+ A) ^ A) 
V > 0(t)-(0(s)+A) t-s 



where we used that for <j>(t) — (4>(s) + A) ^ A, tttt^ h , A) < 73- in the last inequality. Taking i? = R /6 we 



obtain that V(xi,x 2 ) G Bn(s,t,x,x'), d^ t ^_^^ +A ^(xi, x 2 ) ^ 2i? - We therefore derive from (4.8) and (4.10) 
that 3c > 0, 



(M,x,x') ^ Co\<j>(s) + A - s)- d / 2 (t - 0(t))- d / 2 (0(t) - (>(«) + A))- d / 2 / I(s 1 ,s a )eBi l («,t,*,*')<ki<k2. 



/(K d ) 2 

Since - (0(s) + A) ^ i - s and there exists c > s.t. \{x 1 G M d : c^ (s)+A _ s (a:, x i) ^ ^11 ^ 5 (0(s) + A - 
s) d / 2 , |{x 2 G M d : d 2 _^^ (x 2 , x') < J?}| ^ c(t — <f>{t)) d / 2 where |.| stands for the Lebesgue measure of a given set 

in R d , we derive (4.9) from the above equation up to a modification of Co- 

It now remains to do the chaining when for < j < j' ^ N, (x, x') G (R d ) 2 we have d 2 ^_ t . (x, x') > 2i? ^ 1- 
Set L = \Kd 2 . l _ tj (x 7 x')~\, for if ^ 1 to be specified later on and h := (ty — tj)/L. Note that L ^ 1. For all 
i G [0, L\ we denote s, = tj + ih,yi = x + j=(x' — x) so that so — tj, — tj> , y$ — x, %/l — x' . Introduce now 
p := d tj ,- t:i (x,x'){t j i -tjY^/L = \x' -x\/L and for all i G [1,L- 1], := {x G M d : |x-t/i| < p}. Note that 
with the previous definitions Vi G [0, L — 1], \y%+\ — y»| = |x' — x|/Z/ = p. Thus, 

Vxi G -Bi, I a; - xi| < 2p, Vi G [1,L - 2], (xi,x i+ i) G B { x 1^ - a; i+1 | ^ 3p, 

Va;i_i £ B L _i, |a?x,_i - ar'| < 2p. 

We can now choose K large enough s.t. 

3 P /Vh = 3d t .,- tj {x, a/)/VZ ^ (i?o/12) 1/2 (4.12) 

so that according to (4.9), denoting x = x,xl = x' , for all i G [0, L — 1], G R d ,(xi,x i+ i) G Bi x 
p A ' v (si, Sj+i, 2;^, x i+ i) ^ CQ V hT d l 2 (with the slight abuse of notation B = {x},i?L = {x'} andp A ' y (0, /i, x , 
xi) =p A (0, /i,x,xi)). 
We have 

p^{t h t y ,x,x') >E tj)X \l n L lx&e p A ' x ^-0( 3L _ 1 ,t j ,,X^ L _ 1 ,x')\ . (4.13) 
To proceed we have to distinguish two cases: h > A and h < A. 



NON ASYMPTOTIC BOUNDS FOR THE EULER SCHEME 



17 



If ft ^ A, write from (4.13), 

p A (t j ,tf,x,x / ) > E ( , 



I n L-i x ± eB E[p 



A ' X <H'L-l) 



{sL-i,tj> , , x')\X A Li , Xf {sL _ i) ] 



Since we consider the events X A L1 E B L _\, we derive from (4.11), (4.12) that \X A L1 - x'\/Vh ^ 2pj\fh 
3d t .,- tj {x,x')/y/L sC (i? /12) 1/2 . Hence, from (4.9) 



p A {t j ,t jl ,x,x l ) > c^h-^Et 



= c^h~ d l 2 E u 



tj ,x 



Lt-2 vA,-o PfX A , e B L AX A ^X A 



SL-2 ' (j>(sL-2) 



Now P[X S A L1 g B l ^\X a _ 2 ,X a sl 2) } = / Bt _ iP 1 ^-^(s L - 2 ,s L - 1 ,X A _ 2 ,y)dy, but since we restrict to 

X A _ 2 g Bl-2, according to (4.11), we have for all y g |^_ 2 - y\/Vh ^ 3p/Vh (i? /12) 1/2 for the 

previous i? and therefore (4.9) yields 

p A (t j ,t r ,x,x') ^ {^h-WftBL-^^L-ix^g 1. 

Iterating the process we finally get 

L-l 

p A (i j ,V,^,a ; ')^(c - 1 ft- d / 2 ) L ni B il- 

i=l 

Observing that 

35 > 0, Vz g [1, L - 1], m > (4.14) 
we obtain from the previous definition of ft and L: 

p A (t J ,t f ,x,x / ) > (c^h- d / 2 ) L (c P d ) L - 1 (4.15) 
> c - 1 (V-^.)- d / 2 exp((L-l)log(c - 1 c(p/^) d )) >ct 1 {t j ,-t j )- d ' 2 exp{-a% jl _ ti {x,a/)) 

for a suitable c up to a modification of c . 
- If ft < A. We have to introduce for all k g ^ : = e 1°, £- 1], «z € [t fc , [}■ Rewrite from (4.13) 

p A fe,V^,x') ^ E*„,[I n /- ln ^^.P^^isL-utj^X^x')}. 

k=j zei k »( 1 

Define for all fc g [.?,./), i g [1, tt-ffc], I % k g /fe and < S/i < S/2 < • • • < s / «i fc < tk+i- In particular, for all 
i € [1, jj/fe — 1], Sji+i — Sji = ft. Rewrite now, 



p A (tj,t r ,x,x') ^ E t . iX [I n ^_ 2 ^^ x ^ eB E[I nisi ^^ x ± eB .p ' t i'- 1 (a L - 1 ,t r ,Xf l L _ 1 ,x')\J : s tIji _ 2 ]]. 

(4.16) 
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Introducing 



" J 
i'-2 



X 

i=2 



j'-l J ,/ lli=2 B J>, J 

i'-i 

AX A X ir 

|| p V-i( Sji ,Sji+i ,Xi,x i+1 )p ' V-i( s JV-i> Vi^ttV-i' ^ If d^il-^s J. 

i=2 i=2 l'-2 



we derive from (4.11), (4.12) and (4.9) 



^•'-U > (Cb 1 ^ 2 )^'- 1 II 1^/* I / p A '^'- 2 (* ^1)^1 

f A X A 

/ P ' v " 2 (« /V-2' s / 1 ,_ 1 ' X s A »/,_ 
Jb '),_ 1 



Plugging this estimate in (4.16) we obtain 

p A (t 3 ,t f ,X,X / ) > (Co 1 ^-^ 2 )^'-!^)^'-!-! 



xE t . ;X [I ,_ 2 A / p ' V-2( s 2jS x s A ,a:i)da;i] 

<==J »< 'Jb,i V-2 3 Vi 2 

i'-i 



> ( Co - 1 / l - d / 2 )» / ,'- 1 + 1 (cp d )^'-E t ., x [I n , : _ 2n 



using once again (4.12), (4.9) for the last inequality. Iterating this procedure we still obtain (4.15) and can 
conclude as in the previous case. 

Chaining in case (b). If d 2 _ s (x,x') ^ c _1 i? , for c large enough, we derive similarly to case (a) that 

3c > 0, V0 < s < t < T, Vy, p A ' v (s, t, x, x) > ^{t - s)~ d . (4.17) 

Similarly to the previous paragraph we reduce to the case <f>(t) — (</>(s) + A) > A. Then, equation (4.10) still 
holds and for the previous set Br with the current definition of d 2 (., .). From standard computations, we derive 
taking a suitable R that V(xi,X2) € Bn(s,t,x,x'), d 2 ^^ j _^ 8 ^ +A ^{xi,X2) 2i?o- Therefore, 

P A <y(s, t, x, x ') > c^{<p{ S ) + a - s )- d {t - mr d m) - (m + ^)r d I i {xl , X2)eBR{sM )^idx 2 . (4.18) 

J(R d ) 2 

Define now V(u,y) € (0,T] x M d , i? > 0, 

B R {u, y) := {z G K d : - y x - d ' \ 2 /u < i?/7, - - y^uf/u 3 i?/24}. 
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We have that Vz € Br(u, y): 



| M'_ y M'|2 \ z d'+i,d_ y d'+i,d_ v 1 - d +z lld u \ 2 



<(y,z) ■■= 1 7T — L + 6 

j z d' + l,d_ y d> + l,d_ y l,d' u \2 \ z l,d' _ y l,d'\2 

^ o r / ^ it. 

u 6 2u 

Hence Br(c))(s) + A — s, x) x Bn(t — 4>(t), x') C Bn(s,t, x, x') and therefore 3c > 0, \Br(s, t, x, x')\ > c(t — 
4>{t)) d ((t>{s) + A - s) d which plugged into (4.18) yields (4.17). 

It now remains to do the chaining when d 2 . r _ tj (x,x') ^ 2R . The crucial point is to choose a "good" 
path between x and x' . In the non degenerated case it was naturally the straight line between the two points 
(Euclidean geodesic). In our current framework we can relate d%. r _ t . (x, x') to a deterministic control problem. 
Introduce: 

I{t f -tj,x,x') =inf{ / |</?(s)| 2 ds,0(O) = x, (f>(tj> -tj)=x'}, <f> t = Afa + Btpt, (CD) 
Jo 

with A= ( J d ' xd ' ° d ' xd ' ),B= ( l d ' xd ' ), ip G L 2 (\0,t r -ti],R d '). Problem (CD) is a linear determin- 

\ id'xd' U d 'xd' J \ Od'xd' / 

istic controllability problem that has a unique solution reached for 

^ = B*[R(t r - t h s)]*[Q^_ h ]{x' - R(t r - tj,0)x), (4.19) 

where R stands for the resolvent, i.e. V0 ^ t, to < ty — tj,d t R(t,to) = AR(t,t n ), R(t ,t ) — Idxd and 
Qty-tj — J 3 ' 3 R(tj> — tj,s)BB*R(tj> — tj,s)*ds is the Gram matrix, see e.g. Theorem 1.11 Chapter 1 in 

Coron [Cor07l. For (CD) the resolvent writes R(t,t ) = ( . xd ^ d xd ) and therefore the Gram 

\ [t — tojld'xd' Id'xd' / 

matrix of the control problem corresponds to the covariance matrix of the process X t =x + J AX s ds + BW t at 

time tji — tj , that is Q t ., -t, = ( , \i^^ xd ' l/ 5 ' ^j\^/t\ d ' Xd ' Y Hence, explicit computations give: 

\ (tj'-tj) / 2i d'xd' [tj> ~ tj) /6ld>xd> ) 



w r„ i (x') M ' -X 1 ^' „. (x') d ' +hd -X d ' +1 > d -{x'y^'iti, -ti), 

v s e [o, t r - tj], <p s = l^___[ 6s _ 2 ( tj , - tj)] + 6^ (t r -t-y — [tjl - tj - 2s] ' 

(4.20) 

and thus, \lity — tj,x,x') — d 2 t _ t .{x, x') defined in (4.7). Now we have a candidate for a deterministic curve 
around which we can do the chaining. It is simply the deterministic curve (4>s)selo,t jl -t j ] solution of (CD) for 
the above control (<Ps)se[o,t j ,-t j ]- 

To complete the proof of the chaining it remains to specify how to define the (sj)ij>i, (yi)i^i and the associated 
sets. Recall that 2i?o ^ 1- We set here L := \Kd 2 , r _ t .(x, x')~\ for an integer K ^ 3 to be specified later on. 
In term of the new distance, L is similar in its definition to the one of the previous paragraph. Define sq = 0, 

s, :=inf{te [at-utj'-tj] : / |<^ s | 2 ds = I(t r - tj, x, x')/L} A (s;_i + (t r - tj)/L)l ai _ 1<t tj){1 _ 2 } + (tj. - 

ij)I Si _ 1 ^(t.,_t i )(i_ 2.)) i ^ 1- The previous conditions on Ro,K give the well posedness of this definition. 

Lemma 4.2 (Controls on the time step). Set for all i > ; Si := s i+ i — Sj. There exist a constant c\ < 1 a^rf 
an integer L e [L — 1, L/ci], s.i. S£ = tj' — tj and 

Vie [0,1-2], Cl !2l_^^ £jSS !2l_^ 7 t -L_l±^ £L _ 1 ^ 2 t .L_l±. (4.21) 
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Proof. We first set L = inf{fc > 1 : Sk = ty — tj}. The set {k > 1 : Sk = ty — tj} is clearly non- 
empty. The upper bound in (4.21) then follows from the definition of the family (si),j>i. Suppose now that 
Si < (tj> — tj)(l — 2/L) for a given < i < L — 2. Assume also that Sj+i — s, < — tj)/L (otherwise 
e 4 = (^v - tj)/L). Then, J s S4+1 |<^| 2 ds = I(t y - tj,x,x')/L. From (4.19), (4.20), we deduce that 

3c 2 > 0, sup \<p 3 \ < c 2 {tj> -tjY 1 l 2 d t . l -t ] {x,x'). 

0<s<t,/-t,- 3 



Hence, we obtain 



2 I{t r -t j ,x,xf) , dl,_ tj {x,x') 

\ip a \ ds .- z < C2ej ____. 



Recalling that /(tj/ — tj,x,x') = 2d 2 i _ t .(x,x'), the lower bound in (4.21) follows for all i s.t. Sj < T(l — 2/L). 

The bound for L and the last time step are then easily derived. □ 
Define now for all i G [0, L\, yi — (j) Si (in particular y = x and yi = x'), and for all i G [1, Z — 1], 

:= {z G M d : \Q^[R[s u flj-Ow-i - z)| + |Q^ 2 (« - < 2i? ^ 1/2 }, 

where p := dt,-t j {x, x')(tji — tj) 1 ' 2 / ' L. Because of the transport term, we are led to consider sets that involve 
the forward transport from the previous point on the optimal curve and the backward transport of the next 
point in the above definition. Equation (4.13) still holds with L replaced by L. Following the strategy of the 
previous paragraph concerning the conditioning, the end of the proof relies on the following 

Lemma 4.3 (Controls for the chaining). With the previous assumptions and definitions we have that for K 
large enough: 

Vie [1,1-2], \/(xi,x i+1 )eBiXB i+1 , 

dl (xi,x i+ i) sC 2R , 

9 , , (4-22) 
Vsi G B u d 2 Si {x, Xl ) < 2R , 

Vzi_! G Bi^dl^xi^x 1 ) ^ 2R . 

For the same c\ as in Lemma 4-2, 

ViG [1,1-11,1^1 >c lP 2d , (4.23) 

where \Bi\ stands for the Lebesgue measure of the set Bi. 

Indeed, exploiting, (4.21), (4.22) (rcsp. (4.23)) instead of (4.11), (4.12) (resp. (4.14)), the proof remains 
unchanged. The proof of Lemma 4.3 is postponed to Section 4.3. 



4.3. Proof of the technical Lemmas 

4.3.1. Proof of Lemma (4.1). 

The key estimate is the following control of the convolution kernel H A . There exist c > 0, C ^ 1, s.t. for all 

< j < f < N, x, x' G R d , 

\H A (tj,tj,,x,x')\ s$ C(tj> - tj)- 1+v/2 Pc (tj, -tj,x,x'). (4.24) 
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Indeed this bound yields that for all j < f ^ N, x, x' E M. d 

j'-j-i 

\p A ® A H A (tj,tj>,x,x')\ < A ^2 / p A (t j ,t j+k ,x,u)\H A (t :j+k ,t : j l ,u,x')\du 

k=o 

C 2 A ^ (V - t j+k )- 1+ ^ 2 p c {tj' - tj,x,x') < C 2 (t r - t 3 f' 2 B (l, |j Pc (t r - tj,x,x') 

k=0 

using the inequality p A {tj+ k ~ tj,x,u) Cp c (tj ,tj +k , x,u) (cf. Lemma 3.1 of [KMM09] in case (b)) and the 
semigroup property of p c for the last but one inequality. The bound (4.5) then follows from the above control 
and (4.24) by induction. 

Proof of (4.24) We consider two cases. 
- j' = j + 1. From (4.3) we have in this case, Vrr, x' e R d , 

H A {tj , t, y ,x,x') = A- 1 (p A - p A ) {tj , t r ,x,x') 

which are Gaussian densities. In case (a) we have 



r a , , , ., x : ( G((VKa(t j ,x))- 1 (x f -x-b(t j ,x)A)) _ GdVAait^x'jr^x' ^ x - b(t j7 x')A)) 

(A d det(a{tj,x))) 1/2 (A d det(a(^,x'))) 1/2 



H A {tj,t r ,x,x') = A" 



where Vz e R d , G(z) = cxp(— |z| 2 /2)(27r) d l 2 stands for the density of the standard Gaussian vector of R d . In 
case (b) we get 



H A (t j ,t r ,x,x') = A- 1 (2Vi) d ' x 



( (( {AV 2 o-(t j ,x))-\(x>Y> d ' -x^ d ' -h(t j,x)A) 

U 2V3(A 3 / 2 a(t 1 ,x))- 1 ((x / ) d ' +hd - ^ d ' +1 ' d - xl ' d '+f> 1 ' d ' A) 



A d det(a(t,,x)) 



V 



// (A^ait^ixT))- 1 ^) 1 ''' -x 14 ' -b,{t h x')A) \\ \ 

\\ 2V3(A 3 / 2 <j(t 3 ,(x') A ))- 1 ((x') d '+^ d - x d '+^ d - g^jg )^' A ) J j 
A d det(a(i i ,(a;') A )) 



where (x') A := x' — ^ (x') 1 *^ A ) auows to have good continuity properties to equilibrate the singularities 

coming from the difference \x - (x') A \ ^ |(a;') M ' - x^ d '\{\ + f ) + \(x') d ' +1 ' d - x d ' +1 ' d - ^—±^^-A\ with 
the terms appearing in the exponential. In all cases, tedious but elementary computations involving the mean 
value theorem yield that 3c > 0, C ^ 1 s.t. 

\H A (tj,t r ,x,a/)\ ^ CA- 1+r '/ 2 Pc (A,x,x'). 
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- / > 3 + 1- We write in case (a) 

H A (t 0l t r , x, x 1 ) = A ^J G i z ) { (p A fe+i> tj> , x + 6(tj, a;) A + A 1 / 2 a(t j ,x)z, x 1 ) - p A {t j+u t r ,x, x'fj 
- (p A {t j +i,t j/ ,x + b(t j ,x')A + A 1/2 a(tj,x')z,x') -p A (t j+1 ,tj>,x,x'))}dz := 7\ (a) -T 2 (a) . 

Now exploiting that j Rd G(z)zdz = 0, a Taylor expansion at order 3 of T^.T^ yields 

H A (tj,tj>,x,x') = {b{tj,x) - b(tj,x'),D x p A (t j+1 ,tj>,x,x')) + ^Tr^(a(tj,x) - a(tj , x'))Dlp A (t j+1 ,t r , x, x')j 

+R A (tj,t r ,x,x') := (H + R A )(tj,tj',x,x'). (4.25) 

In the above equation H is the difference of the infinitesimal generators at time tj of the processes (X t ) t ^ 
satisfying (1.1) and the Gaussian process X t = x + J t b(s,x')ds + J t a(s,x')dW s , t ^ tj, which can be seen 
as the continuous version of the frozen Markov chain (X A ) ie ^ N ^ introduced in (4.1), applied to the Gaussian 
density p A (tj+i,tji , - ,x') at point x. The remainder term writes 



3A- 1 V f dzG(z) [ d<5(l-<5) 2 

\v\=3 J * d J ° 



R A (tj,tj/ ,x,x') — —Tr^(bb*(tj,x) — bb*(tj 7 x / ))D 2 ,p A (tj + i,tji,x,x') S J + 
D»p A (t j+1 ,tj,,x + S(b(t j ,x)A + a(tj , x)A 1/2 z) , x 1 ) ' x) A + ^ft ' x) 



D v x p A {t j+l , t r , x + 6(b(tj , x')A + a{tj,x')A 1/2 z), x') 



,Ab{t j ,x')A + a{t j ,x')A 1 / 2 z) 1 



using the following notations for multi- indices and powers. For v = (yi, vj) e N d , x = (x\, xj)* set 
\v\ =v x + ... + u d , v\ = v x \...v d \, (x) v = x\\.. x u d d , D v x = D v x \...D v x d d . Recalling the standard control 

3c> 0,C^ 1, \/u,\u\ 4,V0 < j <f < N, {x,x') e (R d ) 2 , \D v x p A {tj,t r ,x,x')\ ^ C{tj,-tj)-^' 2 p c {tj,-tj,x,x') 

(4.26) 

for the derivatives of Gaussian densities, we obtain: 



\R A (tj,tj,,X,x')\ < CA^ity -t ] + 1 )- X p c {ty - tj + \,X,x') + 3A _1 



]T [ dzG(z)[ 1 d5(l-S) 2 



A 



A 



[p v x p A {t j+l ,t jl ,x + 5{b{tj,x')A + (j{t j ,x')A 1 l 2 z),x') - D v x p A {t 3 , t r , x + 5(b(tj, x)A + a(tj, x)A 1/2 z), x'fj x 



{b{tj,x')A + a{t j ,x')A 1 / 2 zy 



Pc(tj> — tj+i, x, x) + 



A 1/2 |x - x'\ v 

(V-t j+1 )3/2 



dzG(z) / ds(i -sy\z\ d x 



Pc (t r ~ t 0+1 ,x + 8{b{t J ,x)A + a{t ,x)A 1 ' 2 z),x') + ^ [ dzG(z) f d6dj(l - 6) 

[tji - tj + i) J R d J[os] 2 



2| z| 4 x 



p c {t r -t j+1 ,x + ^5[(b(tj,x) -b(tj,x'))A+ {a(tj,x) - a{t j ,x'))A 1 l 2 z],x'). 

Now, using the inequality Ve € (0, 1), \x - x' + p\ 2 > \x - x'\ 2 (l - e) + \p\ 2 (l - e~ r ),Vp e M d , taking p = 
5(b(tj,x)A + <7{t ,x)A 1 / 2 z) and p = j5[(b(tj,x) - b(tj,x'))A + {<y{t ,x) - a(tj, x'))A 1 / 2 z] respectively in the 
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first and second integral we get 



ifi a (*,, t f , x, x')\ i _^ )J/2 p ( i-. )c (tj. - t j+u x, x') (i + ^'y jw l/^'if dis<i - ■ s » 2 i= 



Choosing e sufficiently close to one the above integrals are finite and therefore for different c, C depending on 
e as well, we have 

\R A (tj,t f ,x,x')\ ^ C Pc (tj> -t j+1 ,x,x') ^1+ ^ZJ^j ) < C(V -t J+1 )- 1+ "/ 2 Pc (V (4.27) 

Now with the definitions of (4.25) we also have from (4.26) 

\H{tj,tj,,x,7/)\ C{t r -t j+1 )- 1+T <l 2 p c {t r -t j+1 ,x,x'). 
Plugging this last estimate and (4.27) in (4.25) we derive 

\H A {tj,t r ,x,x')\ < C(t r - t J+1 )- 1+n / 2 p c (tj> - t j+1 ,x,x') < C(tj, - tj)- 1+r,/2 p c {t r -tj,x,x'). 
For case (b) we have 

H A (tj,t f ,x,x') = A- 1 / G(z) {(p A (tj +1 ,t r ,x A + B A (tj,x) + Y, A {t 3 ,x)z 1 x') -p A (t j+1 , t r ,x A ,x')) 
- (p A (t J+1 ,t r ,x A + B A {tj,x') + £ A (t jt (i')Ajj') z,x') -p A fe+i, V,x A ,x'))} dz := if 5 - T 2 (b) , 

where s A := x + ^ ^, {x') A ,j,j> ■= x' - ^ ^y,d' ) (*/ - and we set V 2/ e Rd > sA fe,2/) ■= 

( ^,: mm ) ■ -*..» - ( h » )■ - — - - - ■» - 

(a) on Taylor expansions. Actually the d' first component in the above expression can be handled exactly in 
the same way. Anyhow, in order to deal with the d! last components we perform the expansion of , 
around the point xa- We obtain: 

H A {tj,t r ,x,x') = (b 1 (t j ,x)-b 1 (t j ,x'),D xl , d ,p A (t j+1 ,t j >,x A ,x')) (4.28) 
+ l;Tr{(a(t j ,x) -a{tj,(x')A,j,j'))Dl lidl p A (t j+1 ,t r ,x A ,x')} 
+R A (tj,tj> ,xa,x') := (H + R A )(tj, tji, xa, x'), 

where D xl:d > denotes the differentiation w.r.t. the first d! components and similarly to (4.25), H is the difference 
of the generators at time tj of the processes (X t )t^o satisfying (1.1) and the Gaussian process X t = x + 

/** ( (A^) 1 ^' ) ds + ^ Ba ( s ' x ' ~ ( (°')M' ) ~ s )) dH/s ' f G ( continuous version of (A A ) ie[j - n 
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introduced in (4.2)), applied to the Gaussian density p A (tj + i,tj> , -,x') at point xa- The remainder term writes 

R A (t j ,t j ,,x,x') = ^Tr^hbUtj^) - hblitj^'^Dl^^it^tj^XA,^ + 

j (^(^piK{tj^ x ) - hbl(tj,x'))A + (a(tj,x) - a(tj, (x')A,j,j'))^j D 2 xlid , xd , +lid p A (t j+1 , t f ,x A , x')j j 

- 1 / dzG(z) [ d 1 \(D xd , +1 , d p A (t j+1 ,t r ,XA+j(B A (t j ,x) + Z A (t j ,x)z),x'), (B A (t 3 ,x) + X A (t 3 ,x)zf +hd ) 
Jm Jo 1 x ' 

D^^itj+ut'^XA+^B^t^x') + £ A {t h {x')a^)z),x'1 (B A (t J ,x') + £ A {t h (x')a^) zf +1 ' d )}} + 

2A- 1 V [ dzG(z) [ d 7 d ( 5(l- ( 5) 7 2 x 

D D xl , d p A (t j+1 , t, y , xa + 5l(B A (tj ,x) + Z A (tj ,x)z), ^ ((^fe'^+^fe'^)) {B A {tj } x) + S A {t . } x)z) i,d' 



C e ^i, d p A fe+i,V,a;A + (57(B A fe,a:') + S A (^,(a;')A, JJ 0^),a;')- 



0! 



x (B A (i j , a; ') + S A (t j ,(a ; ')A ; ^)^) 



l.d' 



| := (i? 1 ^ + R 2A + R 3A )(t, , t r ,x,x'). 



Let fi = (fi\, ■ ■ ■ , fid') € N d , i/ = (Vi, • • • , Vd') € N d be multi-indices. Similarly to (4.26) we have, 

3c> 0,C > 1, V(/i,i/), |At| < 2,,\u\ 4,V0 < j < ./ iV, (x,x') e K d x K d , 

l^i.*^^^^,*')! < C(V -^.)- (|y|/2+3/2M) Pc(V -*,•,*,*'). 

(4.29) 

The above control yields 

\R 1A {t ,t r ,x,x')\ < Cp c (tj> - t j+1 ,x A ,x') ^ Cp c (t jr -tj,x,x'), (4.30) 
up to a modification of c, C in the last inequality, observing that 

|(*y + M _ x f +M _ C) 1 -' (t ., _ t . +l) |2 _ _ fr'^W* _ tj+i) _ X M' A |2 



(tjv-^+i) 3 fe'-ij+i) 3 

fe' -^+i) 3 

Let us now turn to R 2 ' A (tj,tj>,x,x'). From (4.29) we obtain: 



<S ~ ^P(l-)c(*i' ~t j+ uXA,x') { (jt .,_ tj + l) l / 2 + (t s ,- tj + 1 )$ L dZG{z) ^ X 

exp (c\a\l\z\ 2 (e-i 1)) + ~ ( f A ff' / dzG(z)\z\ 2 exp (c| CT |L kl 2 ^" 1 - 1)) 
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Taking e sufficiently close to 1 the above integrals are finite. Also 

\x - (a/wr < c{\x^ d ' - (z') M 'r + \x d ' +i ' d (x'f +1 ' d (z') M (v - 

s$ C{\x^ d ' - (x'^'r + \x d ' +1 ' d - (x'f +1 ' d - xl " d+ ^ l "\ t r - tj )\" + \x l > d ' - (ar') 1 ' d > (^y^)" ),(4.31) 

and similarly to (4.30), > 1 _^/ 2 P(i-E)cfe' — *j+i> ^A, ar 7 ) ^ Cp c (tj> — tj,x, x') up to a modification of c, C in the 
previous r.h.s. Hence, from the definition of p c in case (b) we deduce: 

|i? 2 > A (^,^v,a;,a:')| ^ C(t r - tj) _1+,?/2 Pc(*j' -ij.ar.a;'). 

Thus, proceeding as in case (a) using (4.31), (4.29) to handle R 3 ' A (tj, t r ,x,x' ), we eventually derive that 
\R A (tj, tj>, x, x')\ ^ City — tj)~ 1+TI / 2 p c (tj, — tj,x,x'). From the definition in equation (4.28) we also obtain 
j-fffo,^/,^')! sC C(tj> - tj)- 1+7 i/ 2 p c (tj> -tj,x,x') using (4.31), (4.29) and the statement (4.24) follows. 

Remark 4.1. Note that the time dependence in the frozen dynamics (4.2) somehow corresponds to the backward 
transport of the terminal condition. It is crucial in order to allow from (4.31) the compensation of the exploding 
terms associated to derivatives in x 1 ^ of order greater than 2 and derivatives in x d +1 ' d of order greater than 
1 appearing in the kernel H A . A similar construction was used in [KMM09]. 

4.4. Proof of Lemma 4.3 

Let us first prove (4.22). We begin with (x i} x i+ i) e Bi x B i+ i, i € [1,Z — 2]. From Section 4.2, one can 
check that \Q7^ 2 {R{si+i, Si)xi - x l+ i)\ 2 = 2d 2 ^x il x i+ i). Hence, 

Qi := d ei (xi,x i+ i) = -^=\Q~ 1/2 (R(s i+1 , Si)xi - x i+1 )\ < c\Q~ 1/2 (xi - R(si, s i+1 )x i+ i)\ 
V2 

< c{\Q7 t 1/2 (xi - R(si, + \Q- 1/2 (R( Sl , s l _i) 2 / l _i -y t )\ + \Q^' 2 ^ - R(si, s l+1 )x l+1 )\} 

■■=Q\ + Q 2 + Ql 

One has 

denoting for all z € := z 1 ^ ,z 2 ■= z d +1,d with a slight abuse of notation. Now, from (4.21), ei/(Kp 2 ) ^ 

ci ((ty -t )/L)/(Kdl^ t] {x,x'){t r -t )/L 2 ) = Cl ^ (rEx<) . Thus, recalling L = fA^ „_ t , (s, a;')! , 3c > 

0, Vje [1,2], ^^cand 

2 

Qj < c^(X 1 / 2 p) 1 - 2 -''|(x J - J R(^,^_ 1 )y i _ 1 ) J -| < clQ-^ixt - ^(ai.si-ijyi-i)! < ci? ^ 1/2 , 

exploiting a;j G Si for the last identity. The term Qf could be handled in a similar way so that Q\ + Qf < 
cRoK- 1 / 2 . Now Q 2 := ^d ei { yi ,y i+1 ) < I( S ,A + i,fc!/ l+ i) 1/2 < c (/ s s ; +1 |^| 2 da) V2 < c ^g^ < 
Hence, for all « e [1,Z — 2], ^ 2i?o for K large enough independent of ty — tj. Eventually, for x\ G 
e B L _ 1 the terms Q := \Q7 1/2 (R(si,0)x - xx)\ and := - tj, s L _ 1 )x L _ 1 - x')\ < 
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c \Qel-^ ( x L-i ~ R( s L-i'tj' ~ tj) x ')\ can be controlled as the previous Q\, i e [1, L — 2] from the definitions of 
B 1 ,B L _ 1 , so that Q t ^ 2R Q , i e {0, L - 1} as well. This proves (4.22). 

It now remains to control the Lebesgue measure of the sets (B i ) ie ^ 1 ^_ 1 j. Define for all i G [1, L — 1], Ei := 

{z£R i : \Q~K^{Vi - 2)1 < 2 J R (3i^ 1/2 ) -1 }- O nc has 3c := c(d) > 0, > cp 2d . Let us now prove E, C B*. 
Write, for all z e Ei, 

Ri ■= \Q^2 2 {R{si,Si-i)yi-i - z)\ + \Q~^^{z ~ R(si,s l+1 )y i+1 )\ 

< |Q^p2 2 (-R(si,Si-i)j/i-i - 3/»)| + 2|Q i ^ / 2 2 (y l - z)| + \Q~^{yi - J R(s J ,s J+ i)y i+ i)| 
:= Rl+Rl + Rl 

The previous definition of ^ gives R\ < 3 ^" 2 . Now, arguments similar to those used to control the above 
{Qi,Qi)ieli,M-2} yield 

Since the term ii| could be handled in the same way we deduce that for K large enough and Ro large enough 
w.r.t. the above c, Ri < 2i? -?^ _1 ^ 2 - Hence Ei C i?i which completes the proof. □ 

References 

[Bas97] R. F. Bass. Diffusions and Elliptic Operators. Springer, 1997. 

[BGL01] Sergey G. Bobkov, Ivan Gcntil, and Michel Ledoux. Hypercontractivity of Hamilton-Jacobi equations. J. Math. Pures 
Appl. (9), 80(7):669-696, 2001. 

[BT96a] V. Bally and D. Talay. The law of the Euler scheme for stochastic differential equations: I. Convergence rate of the 

distribution function. Prob. Th. Rel. Fields, 104-1:43-60, 1996. 
[BT96b] V. Bally and D. Talay. The law of the Euler scheme for stochastic differential equations, II. Convergence rate of the 

density. Monte-Carlo methods and Appl., 2:93-128, 1996. 
[Cor07] J. M. Coron. Control and nonlinearity. Mathematical Surveys and Monographs, 136, AMS, 2007. 

[DM09] F. Delarue and S. Menozzi. Density estimates for a random noise propagating through a chain of differential equations. 

Preprint of the LPMA, 2009. 
[Fri75] A. Friedman. Stochastic differential equations. Chapmann-Hall, 1975. 

[GL08] E. Gobet and C. Labart. Sharp estimates for the convergence of the density of the euler scheme in small time. Elect. 

Comm. in Probab, 13:352-363, 2008. 
[JO09] A. Joulin and Y. Ollivicr. Curvature, concentration, and error estimates for Markov chain Monte Carlo. 

http://arxiv. org/abs/0904 .1312, 2009. 
[KM00] V. Konakov and E. Mammen. Local limit theorems for transition densities of Markov chains converging to diffusions. 

Prob. Th. Rel. Fields, 117:551-587, 2000. 
[KM02] V. Konakov and E. Mammen. Edgcworth type expansions for euler schemes for stochastic differential equations. Monte 

Carlo Methods Appl., 8-3:271-285, 2002. 
[KMM09] V. Konakov, S. Menozzi, and S. Molchanov. Explicit parametrix and local limit theorems for some degenerate diffusion 

processes. To appear in Annates de ITnstitut Henri Poincare, Serie B, 2009. 
[Kol34] A. N. Kolmogorov. Zufalligc Bcwcgungcn (zur Theorie der Brownschen Bewegung). Ann. of Math., 2-35:116-117, 1934. 
[KS87] S. Kusuoka and D. Stroock. Applications of the Malliavin calculus. III. J. Fac. Sci. Univ. Tokyo Sect. IA Math, 34:391- 

442, 1987. 

[Led99] M. Ledoux. Concentration of measure and logarithmic sobolev inequalities. Seminaire de Probabilites XXXIII. LNM, 
1709, page 120-216, 1999. 

[MS67] H. P. McKean and I. M. Singer. Curvature and the eigenvalues of the Laplacian. J. Differential Geometry, 1:43-69, 1967. 
[MT06] F. Malrieu and D. Talay. Concentration inequalities for euler schemes. Monte Carlo and Quasi-Monte Carlo Methods 

2004, H. Niederreiter and D. Talay (Eds.), pages 355-37, 2006. 
[OV00] F. Otto and C. Villani. Generalization of an inequality by Talagrand and links with the logarithmic Sobolev inequality. 

J. Fund. Anal., 173(2) :361-400, 2000. 



NON ASYMPTOTIC BOUNDS FOR THE EULER. SCHEME 



27 



[She91] S. J. Shcu. Some estimates of the transition density of a nondegencrate diffusion Markov process. Ann. Probab., 19- 
2:538-561, 1991. 

[SV79] D.W. Stroock and S.R.S. Varadhan. Multidimensional diffusion processes. Springer- Verlag Berlin Heidelberg New- York, 
1979. 

[TT90] D. Talay and L. Tubaro. Expansion of the global error for numerical schemes solving stochastic differential equations. 
Stock. Anal, and App., 8-4:94-120, 1990. 



